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Abstract 

We present a calculation of the two-loop anomalous dimension for the transversity distri- 
bution hi(x,Q 2 ), 7^ 1 - ) , in the MS scheme of the dimensional regularization. Because of the 
chiral-odd nature, hi does not mix with the gluon distributions, and thus our result is the 
same for the flavor-singlet and nonsinglet distributions. At small n (moment of hi), 7^ 1 ^ is 
significantly larger than r )^ 9 ^ 1 (the anomalous dimension for the nonsinglet fi and gi), but 
approaches r y^ 1 ' very quickly at large n, keeping the relation > 7^ ,9< - 1 ' ) . This feature 
is in parallel to the relation between the one-loop anomalous dimension for f\ (gi) and hi. 
We also show that this difference in the anomalous dimension between hi and gi leads to a 
| drastic difference in the Q 2 -evolution of those distributions in the small x region. 
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The transversity distribution hi(x, Q 2 ) is the third and the final twist-2 quark distribution 
function of the nucleon. It measures the probability in a tranversely polarized nucleon to 
find a quark polarized parallel to the nucleon spin minus the probability to find it oppositely 
polarized. It thus supplies an information about the nucleon spin distribution not provided 
by the gi distribution. Because of its chiral-odd nature, hi can not be measured in the 
totally inclusive deep inelastic scattering (DIS), but it appears as a leading contribution in 
the tranversely polarized nucleon- nucleon Drell-Yan process |], [2|, [3|, [5|, and the semi-inclusive 
c3 ' DIS which detects pions 0, A's 0, |J, and correlated two pions || in the final states. 

Since hi is twist-2, its Q 2 evolution is described by the DGLAP equation f7j. The lead- 
ing order (LO) splitting function for hi has been known for some time0. In the recent 
literature, three papers || || [l(| discussed the next-to-leading order (NLO) Q 2 evolution 
of hi. Vogelsang [g] carried out the light-cone gauge calculation of the two- loop splitting 
function [[13] for hi and converted it into the anomalous dimension. The present authors pre- 



sented the Feynman gauge calculation of the two- loop anomalous dimension ||. The final 
results of these two calculations in different formalisms agreed. These results were subse- 
quently confirmed by [|l(J in the same formalism as |J. These studies completed the whole 
list of the NLO anomalous dimension for the nucleon's all twist-2 distributions together with 
the known result for f x gl], gj, and g x 0. 



In this contribution, we discuss the NLO Q 2 -evolution of hi following our paper M. We 



1 In Proceedings of the Workshop "Deep Inelastic Scattering off Polarized Targets: Theory Meets Exper- 
iment", DESY-Zeuthen, Germany. Sep. 1-5, 1997. (Ed. J. Blumlein et al. DESY 97-200) pages 157-166. 
(presented by Y. Koike) 
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also present a numerical solution of the NLO evolution equation for hi, comparing it with 
that of the gi distribution. This latter part is a new result not presented in ||. 



2 Twist-two operator for hi(x,Q 2 ) and its renormaliza- 
tion 

The transversity distribution hi(x, Q 2 ) of the nucleon is defined by the following relation ||: 

^e lXx (PS\ij(0)a^ l5 ij(\n)\ Q \PS) = 2 [hi(x, Q 2 )(S ±tl p u - S ±vPli )/M 

+h L (x, Q 2 )M(p^n v - p u n^(S ■ n) + h 3 (x, Q 2 )M(S ±ll n v - S ±v n^)\ , (2.1) 

where \PS) is the nucleon (mass M) state specified by the four momentum P and the spin 
vector S, and two light-like vectors p and n are introduced by the relation, P M = p^ + ^|-n M 
and pn = 1. is decomposed as >S M = (S • + (S ■ p)n^ + S±. In (|2.1|), the gauge 
link operator to ensure the gauge invariance is suppressed, hi, and h 3 are twist-3 and -4 
distributions which are not in our interest here. The quark distibutions ^i,l,3 i n ( |2.1| ) are 
defined for each quark flavor. The replacement of ip — > Ctp T and — > — ?/> T C _1 defines 
anti-quark distributions hi(x,Q 2 ), and it is related to hi(x,Q 2 ) by the relation hi(x,Q 2 ) = 
—hi(—x, Q 2 ). Taylor expanding the both sides of ( |2.1|) with respect to A, one can relates the 
n-th moments of hi t L to the following local operator: 

Vw-Mn = Sn'ipi^a^iD^ . . . iD^J), (2.2) 

where S n denotes the symmetrization among u, \ii, ■ ■ ■ , \x n , and the covariant derivative = 
8^ — igT a A° il restores the gauge invariance. In particular, the n-th moment of hi is related 
to a tower of the twist-2 operators: 

MnlhiQ 2 )] = jf 1 dxx n h 1 (x,Q 2 )= £ dxx n [hi(x, Q 2 ) + (-l)" +1 /ii(x, Q 2 )] , (2.3) 

(PS\e^ 1 ... l ^(Q)A^A^--.A^\PS) 

= ^M n [hi{Q 2 )) ^P n+1 - P,SP n + -J—M 2 A^SP n ^ , (2.4) 

where 9^ Ufll ...^ n is defined as the traceless part of M i, M1 ... Mn defined by the condition, 

y u fiu/i!...^ y u ixvnx...nn y u ^u^ 1 ...^ n u ' \^-°) 

and we introduced a null vector A M (A 2 = 0) which symmetrizes the indices u,p,i, ■ ■ • , u n 
and kills off most of the trace terms in 9^^...^ as usual. 

For simplicity we symbolically use the notation O n for 9^^...^ in what follows in this 
section. The bare-(Ojf) and the renormalized- (O n («)) composite operators are related by 
the renormalization constant Z n {p) for O n as 

OM = Z-\M. (2.6) 



2 



The scale dependence of n (p) is obtained by solving the renormalization group equation 

dO n (jj) 



+ in{g{n))O n {n) = o, 



where y n (g(n)) is the anomalous dimension for the operator O n (/i) defined as 

d 



This equation is solved to give 

0„(Q 2 ) = 0„(^ 2 )exp 



g(v 2 ) p(g) 



(2.7) 



(21 



(2.9) 



Up to next-to-leading order, the anomalous dimension 7 n (/x) and the beta function [3(g) can 
be expanded as 

9 -A 1 ) 
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(2.10) 
(2.11) 



where the coefficients of the /3-functions are well known; (3$ = 11 — (2/3) Nf, f3\ = 102 — 
(38/3) iV/, with the number of quark flavor Nf. Inserting these equations into fl2.9|) , one 
obtains the next-to-leading order evolution equation for the n-th moment of hi as 



A^n[^l(Q 2 )] 
i?n(QV) 

a s (Q 2 ) 

47T 



R n (Q 2 ,/j, 2 )M n [h 1 (/j, 2 )] 

i2 n\ 7i 0> /2/3o 



a s (Q 2 
1 



1 + 



a s (Q 2 ) - a s (^) (3 X y£ 



AO)' 



/3 ln(Q 2 /A 2 ) 



1 - 



47T 

/3ilnln(g 2 /A 2 ) 
/5 2 ln(Q 2 /A 2 ) 



(2.12) 
(2.13) 

(2.14) 



The LO evolution equation can be obtained by putting (3\ and y^ equal to zero in (|2.13| ) 
and ( p.14 ). This NLO effect in hi has to be combined with the NLO effect in the hard cross 
section in the parton level to give a prediction for a relevant physical cross section in the 
NLO level. 



3 Two- loop anomalous dimension for hi(x,Q' 



To obtain the renormalization constants for (J,Vfll '"^ r> in (|2.4| ), we calculate the two- loop 
correction to the two-point function for 9 in the Feynmann gauge with off-shell quark lines 
having the momentum p (p 2 < 0). Those diagrams are shown in Fig. 1. Owing to the 
chiral-odd nature of hi, it does not mix with the gluon distribution. Therefore Fermion loop 
diagrams peculiar to the flavor- singlet distribution do not contribute to the renormalization 
of hi. The calculation of the two- loop Feynman diagrams is rather involved, and we refer 
the details of the calculation to our paper M. 



3 





(b) 









(hi 









(m) (n) 







Fig. 1 : Two-loop corrections to the one-particle irreducible two-point Green function which 
imbeds the operator ■ A. 



In the MS scheme, the final result for the two-loop anomalous dimension (|2.10|) for 

hi reads 

r 



In 



where 



S 2 {n + 1)- 251(71 + 1)- 



+CfCg 



e;s c.79 on 

-165 x (n + l)5 2 (n + 1) - yS 2 (n + 1) + ^-^(n + 1) - y 



C2 /T /"* 



4^(71 + 1) ^ ( ^±±) - 5 a (n + 1) - i ^ - 8S(n + 1) 



+ 9'f n + 1> l 5 ?rr? 1 n 1 (i + (-i) w ) 1 1 

l^~J ~~ 2 S2{n + (n + l)(n + 2) + 4 



35 2 (n + l) -55i(n + l) + 



S'fcfn) 



Si. 



?? 



n i 



n f— lV 
5(n) = £^M(j), 



(3.1) 

(3.2) 

(3.3) 
(3.4) 



and Cf = N ^ N 1 , Cg = N c , Tr = Nf/2 with the number of color N c and the number of 
quark- flavor Nf. Equation ( |3.1| ) shows 7^ 1 ) ~ Inn at large n as was the case for (7^ 
for the nonsiglet fi and g\) Jl2j. To compare their behavior, we plotted 7^'-'"^ for Nf = 3, 5 
in Fie. 2. 
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Fig. 2 : Two-loop anomalous dimension 
7 M(i) f or Nf = 3 (circle) and iV> = 5 
(triangle). 



Fig. 3 : 

iV> = 3. 
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From Fig. 2, one sees clearly that at small n 7^ ( - 1 - ) * s significantly larger than 7^^ but 
approaches very quickly to 7^^, keeping the condition 7^) > 7^^. This feature is the 
same as the one-loop anomalous dimensions 7^ < -°* ) which reads 



y(o) 

In 



2C F 1 - 



(n + l)(n + 2) 





(3.5) 
(3.6) 



and hence 7^°) > 7^ for all n. Actually this tendency is even stronger for the two-loop 
case. We note that for n = the anomalous dimension for fi is zero in all orders because 
of the Ward identity for the vector current. On the other hand, hi projects onto the tensor 
operator i/ja^ip for n — 0, for which there is no conservation law. Hence 7o^°' 1 ' ) 7^ 0. From 
the difference between j^ 1 ' and j!^ 1 ' at small n, we expect that the NLO effect leads to larger 
difference in the Q 2 -evolution in the small-x region. We note that the rightmost singularity 
of 7^°) is located at n — —2 while that of 7^ ( - 1 * ) is at n = —1. Therefore, for the case of 
the transversity distribution, the NLO evolution renders the DGLAP asymptotics at x — > 
compatible with the Regge asymptotics [[15] . 

Since the relevant quantities for the (^-evolution of the moments are /2(3\ — , y^/2f3 
(see Q2.13D ), we plotted it in Fig. 3 for Nf = 3 and 5 cases. From this figure, one expects 
that at small n the NLO effect in the Q 2 evolution is quite different between M. n [hi{Q 2 )} 
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and MnlhiQ 2 )]. Since 7^ ' 1 ) -> as n -> 0, 7^ 1 >/2/3 1 - 7;{ (0) /2/? abruptly drops to zero as 
n — > 0. But this is not the case for h\. It shows a characteristic behavior as shown in Fig. 3. 

As an example of the (Revolution, we show in Fig. 4 (a) and (b) the (Revolution of the 
tensor charge and the first moments, respectively, with the parameters Nf = 3 and A = 0.232 
GeV in ( j2.14 ). (Here we are interested in the NLO effect in the anomalous dimension and 
the /^-function, and thus we adopted the same value for the A-parameter in the LO and NLO 
evolution.) 
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Fig. 4: (a) The LO and NLO (Revolution of the tensor charge normalized at Q 2 = 1 GeV 2 . 
(b) The LO and NLO (Revolution of the first moment of hi(x, Q 2 ) and f\(x, Q 2 ) normalized 
at Q 2 = 1 GeV 2 . 



They are normalized at Q 2 = 1 GeV 2 . At n — 0, only diagrams (e), (f), (i), (j), (k) in 
Fig. 1 survive. They give the anomalous dimension for the tensor charge as 



h(i) 
7o 



-19C 2 



257 52 
9 9 



724 
IT 



104 



N f . 



(3.7) 



From Fig. 4(a), we can compare the LO and the NLO (Revolution of the tensor charge of 
the nucleon. One sees that the NLO effect is sizable as is expected from Fig. 3. In Fig. 4(b), 
we plotted the LO and NLO (Revolution both for Mi[fi(Q 2 )] (nonsinglet) and Mi[h 1 (Q 2 )}. 
Although the NLO effect in the anomalous dimension (the second factor in the right hand 



side of ( 2.13|) ) makes this ratio smaller, the NLO effect in the coupling constant (the first 
factor in the right hand side of (|2.13|) ) completely cancels this effect. For fx, the latter effect 
is actually larger than the former effect. As is shown in these examples, the NLO effect on 
each moment is not large. We will see, however, that it is significant for the distribution 
function itself in the next section. 
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4 Q 2 -evolution of h\ 



To see a generic feature of the actual NLO Q 2 evolution of hi, we have applied the result 
of the two- loop anomalous dimension to a reference distribution function of h\. In the 
nonrelativistic kinematics, h\ = gi, although the motion of quarks in the nucleon is not 
expected to be nonrelativistic. A bag model calculation suggests hi is not very different 
from gi at a low energy hadronic scale. We therefore assume that hi(x,Q 2 ) is identical to 
a gi(x, Q 2 ) distribution of Gluck, Reya, Stratmann and Vogelsang (GRSV) [|l(| (standard 
scenario) at a low energy scale. 

We first recall the general formalism of the NLO Q 2 -evolution. Corresponding to the 
splitting of a quark into an anti-quark, there are terms proportional to (— l) n+1 such as 
those caused from S' k and 5* in the anomalous dimension (|3.1| ). One can thus write for 
in (fTIJ) as 



Rn{Q 



V) 



K q (Q 2 ^ 2 ) + (-i) n+1 Rf(Q 2 ,» 2 ) 



(4.1) 



We solve the evolution equation ( 2.12|) with (2J3) and ( [4.1|) for n = even and odd cases by 
analitically continuing them to all n. This gives the Q 2 -evolution for the two combinations 
hi(x,Q 2 ) ± hi(x,Q 2 ). Numerical solution was obtained following the method of We 
have checked that the method works sufficiantly accurately in the x > 0.05 region for the 
NLO evolution. 
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Fig. 5: The LO and NLO Q 2 -evolution for the w-quark transversity distribution 5u, starting 
from the same input function. 
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Fig. 6: (a) Comparison of the LO (Revolution of Su(x, Q 2 ) and Au(x, Q 2 )- (b) Comparison 
of the NLO (Revolution of 8u(x, Q 2 ) and Au(x, Q 2 ). 
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Fig. 7. Comparison of the NLO (Revolution of Su(x, Q 2 ), Sd(x,Q 2 ) and Au(x,Q 2 ). 

To see the magnitude of the NLO effect, we first applied the LO and NLO evolution of 
hi to the same input function. As an input function we have chosen here the GRSV LO g\ 
distribution function at [x 2 = 0.23 GeV 2 , which is the starting distribution of fl6[| . Figure 
5 shows the result for the u-quark transversity distribution Su(x, Q 2 ) evolved to Q 2 = 20 
GeV 2 . The NLO effect is significant especially in the small-x region. The same tendency 
was observed for the other transversity distributions. 

To compare the (Revolution of the transversity distribution and the gi-distribution, we 
show in Fig. 6 (a) and (b), respectively, the LO and NLO (Revolution of 5u and Am (g 1 for 
u-quark), starting from the GRSV LO and NLO input function for Au given at Q 2 = 0.23 
and 0.34 GeV 2 . Although the LO evolution leads to a significant difference between Su 



and Am at Q 2 = 20 GeV 2 , this feature is even more conspicuous in the NLO evolution, 
in particular, in the small-x region. This is expected from the behavior of the anomalous 
dimension. (Although the evolution of Am shown in Fig. 6 is receiving some effect from the 
mixing with the gluon distribution, flavor-nonsinglet and valence distributions also show the 
same tendency.) 

In Fig. 7, we compare the NLO Q 2 -evolution of 8u, 8d and Am, starting from the same 
input distribution function (GRSV NLO input function for the sea quark distribution for 
gi). The difference between 5u and Am is again significant. Although the input sea quark 
distribution is assumed to be flavor symmetric (5u = 5d at Q 2 = 0.34 GeV 2 ), the NLO 
evolution leads to a violation of this symmetry because of the presence of R^(Q 2 ,fi 2 ) in 
(|2.13| ). However, this effect is very small as can be seen from Fig. 7. | 



5 Summary and outlook 

In this work, we have carried out a Feynmann gauge calculation of the two-loop anomalous 
dimension for the transversity distribution hi(x,Q 2 ) in the MS scheme of the dimensional 
regularization. This completes the calculation of the anomalous dimensions for all the twist- 
2 distributions of the nucleon in the NLO level. We found ^ is significantly larger than 
the j^' 9 (nonsinglet anomalous dimension for fi and g\) at small n, but approaches very 
quickly to 'jl' 9 ^ at large n, keeping the condition > 'jl' 9 ^. The NLO effect in the 
anomalous dimension leads to a drastic difference in the (5 2 -evolution between hi and g\. 
The NLO parton distribution function has to be combined with the NLO short distance cross 
section in the same factorization scheme to give a prediction for a physical quantity. For the 
transversely polarized Drell-Yan process, the latter has been calculated in H . Combining 



these results, the NLO calculation of the double spin asymmetry for the Drell-Yan process 



has been presented in a recent work [18]. We hope the peculiar feature of the hi distribution 



discussed in these works will be measured in the future collider experiments. 
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